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ABSTRACT 


The  purpose  of  this  report  it  to  analyse  the  stiffness-  tod  damping  pro¬ 
ps  rtie*  p£  fluid  film  Journal  bearing*  and  to  determine  the  force  transmitted 
to  the  bearing  support.  The  results  are  presented  in  practical  design  charts 
and  are  given  in  dimension! ess  form  to  make  them  applicable  to  wide  ranges 
of  geometrically  and  dynamically  similar  units. 

The  analysis  assumes  that  the  rotor  vibrations  are  of  small  amplitude, 
Thereby  the  non-linear  oil  film  force  is  replaced  by  gradients,  denoted 
spring  and  damping  coefficients.  The  numerical  values  of  these  coefficients 
are  obtained  by  computer  calculations.  Results  are  given  for  3  bearing  types! 
the  plain  cylindrical,  the  4-axiftl  groove  and  the  elliptical  bearing.  Using 
these  bearings  to  support  a  symmetrical  two-bearing  rotor  the  force  trans¬ 
mitted  to  the  bearing  pedestals,  due  to  a  rotor  imbalance,  is  calculated. 

Thus  the  bearings  can  be  compared  for  a  given  rotor  and  for  known  operating 
conditions  and  the  bearing  with  the  optimum  force  attenuation  can  bo  selectedv- 
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INTRODUCTION 


The  hydrodynamic  oil  film  force  is  obtained  from  Reynold**  equation, 
assuming  constant  viscosity.  The  equation  is  approximated  by  a  finite 
difference  equation  and  solved  on  a  computer.  The  resulting  oil  film  force 
is  a  non-linear  function  of  the  eccentricity,  the  attitude  angle  and  the 
corresponding  velocity  components.  The  non-linearity  implies  a  complex 
relationship  between  the  rotor  and  its  bearings  such  that  in  an  exact  analysis 
the  rotor  and  the  bearing  cannot  be  studied  separately,  but  must  be  analysed  as 
a  system.  Even  if  an  exact  solution  was  available  it  would  not  be  too  useful 
for  design  purposes,  firstly  because  of  the  vast  number  of  bearing  and 
rotor  parameters,  and  secondly  because  it  would  be  almost  hopeless  to  tie 
the  results  in  with  the  supporting  structure.  In  the  present  analysis,.  therefore, 
the  oil  film  force  is  linearised  by  replacing  it  with  it#  gradients,  mathematically 
expressed  in  the  first  order  Taylor  expansion.  This  is  a  justified  approxima¬ 
tion  when  it  is  assumed  that  the  journal  motion  is  email.  A  linear  bearing 
force  vastly  simplifies  the  rotor  analysis  and  makes  it  possible  to  assign 
an  impedance  to  the  bearing,  a  necessary  presupposition  for  any  overall 
investigation  of  the  rotor  and  its  supporting  structure. 

Three  bearing  types  are  studied:  the  plain  cylindrical,  the  4-axial 
groove  and  the  elliptical  bearing.  The  configurations  are  shown  in  figure  1, 

2  and  3.  The  oil  film  force  gradients  are  calculated  on  the  computer  as 
shown  in  table  1-4  and  introducing  the  numerical  values  into  the  linearised 
expression  for  the  oil  film  force  (see  eq.  (6)  and  (7),  page  4  )  yields  the 
bearing  spring  and  damping  coefficients  as  shown  in  table  5. 

Although  the  thus  obtained  data  are  completely  sufficient  for  a  rotor 
vibration  calculation  they  are  not  in  a  too  convenient  form.  The  reason  is 
that  8  coefficients  are  obtained  whereas  the  normal  rotor  calculation  is  set 
up  for  only  4  coefficients,  a  spring  and  damping  coefficient  in  two  mutually 
perpendicular  directions.  No  provision  is  made  for  taking  into  account  the 
additional  4  cross -coupling  coefficients,  Therefore,  it  is  a  matter  of 
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practical  importance  to  eliminate  them.  Unfortunately,  they  are  unsymmetrical 
and  do  not  vanish  by  the  introduction  of  principal  axil.  Instead  another  method 
i»  employed  making  use  of  the  fact  that  linear  bearing  forces  result  in  har¬ 
monic  rotor  motion.  Thus,  it  is  possible  to  replace  the  original  8  coefficients 
with  4  equivalent  coefficients  that  will  give  exactly  the  same  rotor  motion. 

It  it  clear  that  such  a  reduction  depends  on  the  rotor.  A  symmetrical,  two- 
bearing  rotor  is  selected  since  it  represents  the  most  commonly  used  rotor 
design,  As  a  further  simplification  the  rotor  is  given  only  one  degree  of 
freedom  by  concentrating  the  rotor  mass  at  midspan,  The  simplified  system 
it  shewn  in  figure  6, 

The  procedure  is  as  follows!  a  force -balance  in  the  vertical  and  the 
horizontal  direction  is  set  up  for  the  rotor  combining  the  rotor  inertia,  the 
unbalance  force  and  the  bearing  force  represented  by  the  computed  &  spring 
and  damping  coefficients.  This  results  in  4  equations  in  the  unknown  ampli¬ 
tudes  set  up  in  a  matrix.  By  eliminating  the  8  terms,  representing  the 
cross -coupling  effect,  the  matrix  is  reduced  to  the  same  form  as  the  matrix 
for  a  rotor  with  only  4  spring  and  damping  coefficients  and  no  cross -coupling 
coefficients,  Therefore  the  remaining  terms  in  the  reduced  matrix  are  the 
desired  equivalent  spring  and  damping  coefficients.  In  addition,  the  matrix 
is  solved  for  the  rotor  amplitudes,  and  combining  the  amplitude  with  the 
spring  and  damping  coefficients  give*  the  force  transmitted  by  the  bearing. 

All  the  results  are  given  in  dimensionless  form  as  a  function  of  a  dimensionless 
parameter  '#*  •  In  plotting  the  results  "*  is  replaced 

by  the  speed  ratio  m,  see  eq,  (29),  page  14  .  Curves  are  given  for  the 
dimensionless  spring  and  damping  coefficients,  and  the  dimensionless  trans¬ 
mitted  force  as  shown  in  figures  12-43,  The  force  attenuation,  expressed 
as  the  ratio  between  the  aetuair  ransmitted  force  and  the  force  transmitted 
by  a  rotor  in  rigid  bearings,  is  easily  found  from  the  graphs  since  the 
dimensionless  rigid  bearing  transmitted  force  is  also  plotted, 
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THEORETICAL  ANALYSIS 


For  an  incompressible  fluid  Reynold's  equation  may  be  written: 

111  ’  <>**«$  +  Itjffcosd  +  I2e$rtine 

-  6Rco(l-2|)$  +  l2l|f  cose 
where  C+ecose«C  +  €cos(l) 


Introducing:  fi»2Ch  X«2RX  2«Lz  (5*C6 

P-/uN(l-2^®'  (constant  yU  ) 

i  »•  iti 

Reynold*  ^.equation  reduces  to  the  following  dimensionless  equation; 

^^  +  l27rTf4r< 


(2) 


•cose 


The  resulting  oil  film  force  is  then: 

(3) 


K)  U.Ab-U.Agj,  AAAftaa  All  A*  t 

5»/AN(hHi)(ffDd  \  Pcw(a*)«lx<fcc 


whe  re 


(4) 


A-'#1® 


For  the  subsequent  analysis  it  is  necessary  to  linearize  the  force  with  respect 
to  displacement  and  velocity.  The  first  order  approximation  of  the  Taylor 
expansion  will  be  used: 


<*5  -  Je  +  &  id*  * 


W  riting : 


d4 


if*  i  4 

Wi-ii)  ‘  (IW 


and  taking  as  the  reference  for  the  Taylor  expansion  the  steady  state 
equilibrium  position  where 


we  get: 


(5)  a  £  a  Aw[||  dc  +  ^  ed*  +  !|r  d&)  - &d<*] 

To  change  from  polar  to  rectangular  coordinate*  (see  figure  4)i 

XBC?£M#  y»C£S}ttei 

de  «i[eost!  dx  +sln<*dy]  ed«*iE“sM  dx  teosady] 

Since  the  point  )  is  the  steady  state  position*  dx  and  *epra»airt* 

the  dynamic  displacements  and  dx  a  (dx )  and  d|  •-  (dy  }- 

Thus  we  obtain: 

(6)  d!*£M»|(4f^  |a*  $’»»*)  dx  cm  +  ^sftw^di 

+  ($£■  tin* + cm)  dy  +(|||  tin*  dy  } 

on  dfjp  —  Ki,,  dx  - C«„  dx  *  K*,  dy  +  d*y 

i't)  d ^  dx  dx  -  Kw  dy  - d*y 

Three  bearing  configurations  will  be  analysed.  For  this  purpose  equation  (6) 
its  not  in  a  convenient  form  and  must  be  rewritten  in  terms  of  the  force 
components  normally  used  in  bearing  calculations. 

Cylindrical  Bearing 

The  force  is  given  in  terms  of  a  radial  component  F^,  positive  in  the  negative 
radial  direction,  and  a  tangential  component  Ft,  positive  in  the  positive 
A  “direction,  such  that 

ft=~$£oo5*  fjj8 -Fjlsina  +  ftcos* 

c(i- ~  -dfccosd  t  i^sinei  tk  “dFt  sinrt“FtC05d  d& 

d^a  ”dl>5{hd  "“i>cosddd  ■+  d cosoi  “F*$iha  d* 


**4 


I*- 


Expressing  d.F_  and  dF^  by  equation  (6),,  the  coefficients  in,  equation  (?) 

X  U 

can  be  written,!  * 

cos^  +(“£  +  ||) eos<* stool] 

i 

i 

to  C***  &  Xto  [§|j  cos*al  +  ^  slrftfc + -( ^  ^ )  com  s  tod  3 

Kq»iX»  F’l  cosol  -H  sMt  +  (■£  -  II )  tosd  Sim] 

toQa^Aw  [^cos^-JI  §tof<*+  (^*"*  J|t  )cosdSihd] 

(a) 

K^IAm  [H  cos'd  +  stold+  “*  i|“  )  cow- stool  j 

to  i  A«  iJI  cos'd  -  $fh*A+  )  cosd  Stod  ] 

Kif  £Aw  [|  cos’d  + 1|  stoV,  I "l*  H  )  co«d  Si  ltd  ] 

wC***£Aw  [^  cos\*+H  sto'd”  (%  )casds?m] 

4 "Axial  Groove  Bearing 

The  force  is  given  in  terms  of  a  vertical  component  Fv>  positive  in  the  negative 
x-direction,  and  a  horizontal  component  F^,  positive  in  the  positive  y-direction, 
such  that 

f^«“Fv  Fh 

Using  equation  (6)  directly  the  coefficients  in  cq.  (7)  become! 

Ksf  i  Aw  [II  cm  ~  ^  Stoot] 
w  Cxx»  i  Aw  [$fj  eosd  +  ^  stool  ] 
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(9) 


w  Qi,  *  i  Aco  [~||  4  ] 

K^cAa>[^-cos*  -  It l  slh*] 

cosu  4  ^  s'im] 

&>C W*£A«  [~|j|  Sifld4^C0$d] 

Elliptical  Staying. 

The  elliptical  bearing  ia  made  up  of  two  partial  arc  bearings  sailed  the  lower 
lobe*  identified  by  subscript  1,  and  the  upper  lobe,  identified  by  subscript  2, 
The  radial  bearing  clearance  is  taken  as  the  difference  between  the  lobe 
radius  and  the  journal  radius.  Ths  origin  of  the  x,  y~ coordinate  system  is 
located  at  the  bearing  center,  midway  between  the  lobe  centers,  with  the  x- 
axis  vertical  downwards. 

From  figure  5: 

$ae*+to**2£mcos«  e*»  e^hMemcos* 

Sfcbtti* 

Furthermore 

c(e,=  ilcosd.dtf  +5tod(  sty] 

(ll) 

del, «  £  ["Slhtfjdx  4  cysd.dy] 
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defs  4  £*coSflitdx  +  sin**  dyl 


l*dd* «  4  [  Sifted*  +fcdSd®  d)|] 

The  two  lebos  are  calculated  separately  resulting  in  a  vertical  and  a 
hoj?i*ontal  force  component  for  each  lobe,  Them 


These  equations  are  analogous  to  the  4-axial  groove  bearing.  Therefore 
aquations  (9)  are  also  applicable  to  the  elliptical  bearing.  However,  a 
difficulty  is  encountered  in  the  calculation  of  the  derivatives  with  respect  to 
velocity  because  a  pure  radial  velocity  £,  gives  rise  to  both  a  radial  and  a 
tangential  velocity  component  for  the  lobes.  Thus  «~2f  does  not 

hold  for  the  elliptical  bearing  as  it  did  for  the  cylindrical  and  the  4-axial 
groove  bearing,  but  it  is  still  valid  for  each  lobe  taken  by  itself.  Using 
eq.  (5)  and  (?)  together  with  eq.  (12)  we  get! 

Cyr dx  -Cfl d'ij  «  Aw  [  $$j- i  4  -  dle  -  f£j  e,d«,  +  f £  A ] 

Using  equations  (11)  we  get! 

d)CH*«4Xw[||jiCQSdt 

to  C  Xto  [~^j  sihd,  +  ^  sH  +  ^cosd,  -  C0Sdt  j 
toCpa  c  Xfio  [J^  u>H  +  coMt  +  ■+  •^■sihd^j 
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UiC^f  i  Acj[“  sth«,  +  4^  S(hd8  +  ^  COM,  -  *&*•  cose*,] 


Thus  we  may  use  equations  (9)  by  setting: 


(13 )  |||  COSf^t)  +  C05 (oft (ij)  “  ^  Slftfa-4,)  +  ^  S(h  (<*+«») 

,+,^A  ^■CosfdteJ,) 

^,C05(«* "*«*,)  +1^1  cos (oiteij)  —  ^ sln(di-^  f 

Slr»C<A“*<A|)  f-|^  $tb(d+d,)  +  ^C0$(d-el,)  **  ■^eosfo+dj) 

The  above  forces  and  derivatives  o,re  calculated  by  means  of  a  computer  as 
summarized  in  tables  1-3.  The  resulting  spring  and  damping  coefficients  as 
calculated  from  equations  (8),  (9)  and  (13)  are  shown  in  table  4.  These 
results  can  be  used  directly  when  calculating  the  vibrations  of  the  rotor. 
However,  the  usual  calculation  procedure  allows  for  only  4  coefficients,  one 
spring  and  damping  coefficient  in  the  vertical  and  in  horizontal  direction.  No 
provision  is  made  for  taking  into  account  the  4  cross -coupling  terms  Kxy, 
Cxy,  Kyx  and  Cyx*  Therefore  it  becomes  important  to  eliminate  them  to 
reduce  the  original  8  coefficients  to  4  equivalent  coefficients.  Duo  to  the 
non-symmetry  of  the  cross -coupling  terms  they  do  not  disappear  by  the 
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introduction  of  principal  axis'*  Instead,  the  8  coefficients  maybe  combined 
to  4  by  coupling  the  bearing  with  the  rotor  in  such  a  way  that  the  resulting, 
motion  remains  the  same.  This  is  the  purpose  of  the  following  analysis. 

The  rotor  is  a  simple,  symmetrical,  one  -degree -of -freedom  rotor,  It  is 
supported  in  two  identical  bearings  and  the  rotor  mass  is  considered  con¬ 
centrated  at  midspan,  (see  figure  6). 

JLet  0  be  the  steady  state  position  of  the  journal  center  (  i, e. ,  at  aero 
unbalance),  A  is  the  actual  journal  center,  B  ia  the  shaft  center  at  midspan 
and  G  is  the  center  of  gravity  of  the  rotor.  A  force  balance  gives: 

M  4 + k  (xk-  O  *  M  e  <*)x  m  wt 

(w)  *  Mew*6tnwt 

kfyiry*) *  4*- 

The  following  parameters  are  introduced: 


(IS) 

4*fr 

(16) 

Furthermore  the  solution  is  taken  in  the  form: 

Xh*  Aco$cot  +  Bstowt 
h*  ^§~Sr  coswt  +  sihtut 

4a*  E  cos  cut  ■+  Fsthwt 

4*»R  wp^S* cosw^  " 


"9- 


(17) 


'Substituting  eq.  (15)>  (16)  and  (1?)  into  cq!  (14)  yields! 


A 

B 

E 

F 

w 

'■rr 

¥e~ 

xe“ 

(K^x) 

£hC«K 

“Kx, 

1 

~wQk 

(K^HI 

hi  Cxy 

"Kx, 

0 

“■foC^ 

<oG^ 

0 

coC^ g 

“toGjij 

1 

It  is  desired  to  reduce  this  matrix  to  the  same  form  as  a  matrix  for  a  rotor 
without  cross-coupling  terms,  Such  a  rotor  has  only  4  spring  and  damping 
coefficients  which  are  denoted  K^,  Bx,  Ky  and  By,  The  reduced  matrix  is 
then! 


A' 

itwnwwmrw 

J*e 

B 

*€ 

SJ 

X6 

F 

«e 

(K-*) 

CdBx 

0 

0 

1 

( K*-*) 

0 

0 

0 

0 

0 

(K,-*) 

to  Bij 

0 

0 

0 

~t*>  Bu 

(K»-*) 

1 

~  iO  - 


After  ft  substantial  amount  of  algebra,  eq,  (18)  is  reduced  to  eq.  (19)  with 
the  results* 


where* 

“twCy^)  +0(xij'“wC^) 

T||j«  (  Kxk  m  X  •~wC^)  ■+•  ( kiftf+toCsof) 

<2l)  6  *  (k**  -  ac  -  wC^)(K^“««C^)  ■+(k^”'&e  +cuCw^)(K,x'»”^Qx) 

^"(kxx  -*  -U>Ct}s)(k^-'X  +  wC^)  -(k^-wC^lik^+toCim) 


Since  eq,  (19)  are  linear  they  may  also  be  solved  for  the  amplitudes i 

*  W*)'+W 

B  .  *<wM 

*  lK,~i7+  &"BT 

M  eyemw—  ewww>«we>Ma»iiwaweww»iawwiww»e 

*  (k,-*)V((oB,7 

F  ^(K"«) 

*****  (B  »wwwww»ew>  iWww'iiiiiMWi  w  */*■>»»  nrn«itt»ae*ew* 

*  (K,-*mwBi 


(22) 
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The  js-  and  y-amplitudes  are  found  by  substituting  eq.  (22)  into  eq.  (17) s 

C06  (cat1 “fx) 


(23) 


*"**  Of,-*) 


'ij  "v;  t\w  ) 

6jBh 

- -*mn9nm 

Y 


x^  and  y^  may  be  found  similarly  from  ©q.  (17), 

The  force  transmitted  to  the  bearing  pedestal  is  given  bys 

ft »  K/Xe.  +  B/x* 

(24) 

ft*  K^+B^lJ* 

Substituting  eq.  (23)  into  eq,  (24)  yields! 

^■4M  '"to-***) 


’Uny 


r'“  ^ . K 


(25) 


In  the  calculation  the  above  equations  are  mad©  dimensionless  by  dividing 
through  by  in  order  to  make  them  general.  As  an  assistance  in 

plotting  curves  of  the  derived  equations  the  following  auxiliary  expressions  are 
set  up! 
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Fox’  ©© 


Kx~*  K**~  ufQH 
wBx^wCsx+Ksjj 
Ksj™^  HV^wCV 

<26) 

|jjX  *mp&  | 

^  <Hfcs»  J 

I  ”** 

§  “*■ 

Instead  of  expressing  the  rotor  amplitude  in  x  and  y-coordinates  a  better 
physical  picture  is  obtained  by  finding  the  corresponding  elliptical  path  of 
the  Journal  center.  Combining  the  first  and  the  third  of  eq.  (IT)  we  get! 

+« @M4n 

w  +(f)' +(lf]  ■  ^  |[(4),+!f)'  +(ff+(l'Fr-fffi©“®®i 

tqh&A 

Where  a  is  the  major  axis  of  the  ellipse,  b  is  the-  minor  axis  and  ft  is  the 

A  B 

angle  between  the  x-axis  and  the  major  axis,  see  figure  7,  (  -g~),  (*--),  ( 

tp 

and  ( “g~  j  are  given  by  eq,  (22).  Rotor  resonance  may  be  defined  as  the 
speed  where  the  major  axis  is  a  maximum,  This  maximum  is  found  by 
plotting  the  major  axis  as  a  function  of  ft  .  The  results  are  shown  in 
figure  46-47,  From  these  graphs  the  rotor  critical  speed  can  be  found 
directly  for  a  given  rotor  by  a  trial  and  error  process, 
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Thus  the  results  are  obtained  as  a  function  of  ,  but  to  facilitate  the 
interpr station  of  the  results,  &  is  replaced  by  a  speed  parameter.  From 
eq.  (16)  K  in  dimensionless  form  is; 


■  v&U  is  used  instead  of  Ilf  to  present  the  dimensionless  results  as  shown 
in  figures  11-43 ,  When  a  rotor  with  a  dimensionless  stiffness  different  from 
5  is  investigatedi  eq,  (28)  should  be  substituted  into  eq,  (29)  to  find  the  value 
of  (t$j  corresponding  to  the  desired  value  of  .  This  relationship  is 
shown  for  a  wide  range  sf  dimensionless  rotor  stiffnesses  in  fig.  11. 

The  force  attenuation  may  bo  expressed  as  the  ratio  between  the  actual  trans¬ 
mitted  force  and  the  force  transmitted  with  rigid  bearings.  The  following 
relationship  exists: 


SUMMARY  OF  RESULTS 


Effective  spying  coefficient  In  vertical  direction! 

K/  K«  "  ^  ^  +  ^wCn,)] 

Effective  damping  coefficient  in  vertical  digestion: 

w  B* ■  <oQ*+' 4  C  ^  Kxt#  -  CcAiCxy)3 

Effective  so rina  coefficient  in  horizontal  direction; 

Kj*  Kij,  ^  [(J  K(*x  “”^(foC,jr)] 

Effective  damping  coefficient  in  horizontal  direction;; 

w  ■  toCy, "  yji nj  Ify  +  <J 

where ; 

2  J, 

-  x  •*  fuC^)  +  (  Kk,  -^c^) 

(  kj^  ”*3f  CuCjx)  +■  (  kijr  "^fuGn) 

(J  ®  (  kix  if  ki^j  "^c  'I'^jCkh') 

fr)a  ( K**"* - ^C,,)( ^“"if+wCxi,)  -(fcy-tomd^ttoU) 

X*  i  k 

KtKt  ^Qrx  i  Koj  , Co C»«j  ,  kljx  ,  W  C^k  ,  Ktjij  ,  and  are  given  in 

table  5, 
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Transmitted  force  in  vertical  direction 


Px a  mlrt-ft*#) 


Wy« 


toB* 


K» 


Transmitted  force  in  horizontal  direction 


r,  ■  «"|f3§ 


,.  ,.  stnfast- 


wB, 
"  KfS 


Amplitude  of  journal  center  in  vertical  direction 

**”  wAw 

Amplitude  of  journal  center  in  horizontal  direction 


fvfcr  sth("t'^) 


Elliptical  path  of  journal  center 

ilix+&+  E"+  F*J '+ 1  \/w^iv"E*+'F'l‘-'[AF-BEi< 


Major  axis!  a 


Minor  axis:  b  »  V £ W+  Bt+  £*■+  F* J * i  wt+tf+f+Ff- Fa f ~ 

Angle  between  vertical  and  major  axis!  &  «~r“  tan*1! 

2  L(A+BM:*'F* 


•  16. 


whe  re; 


e  is  rotor  mass  eccentricity,  i.e, ,  e 


A 


B 

jffi 

F 


* 


SS 


(K.-kWwBJ' 


eit(Kr^) 

li^*W 


Unbalance  (lb a » in) 
Rotor  Weight  lISsT 
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DESIGN  INFORMATION 


This  section  has  been  prepared  to  assist  the  designer.  The  method  ia 
directly  applicable  and  does  not  assume  familiarity  with  the  underlying 
analysis.  Thus,  the  designer  should  be  able  to  extract  any  desired  infer- 
motion  without  further  reference. 


In  selecting  the  bearing  on  the  basis  of  minimum  transmitted  force,  the 
designer  is  faced  with  the  problem  of  how  to  choose  bearing  type,  unit 
bearing  load,  -ratio,  clearance  and  oil  viscosity.  To  answer  these 
questions,  the  procedure  below  can  be  employed. 

It  is  assumed  that  the  rotor  is  given  such  that  the  following  quantities 
are  known: 

y  &  cl  y 

1)  The  rotor  critical  speed,  COc  /  sec,,  for  rigid  rotor  supports  (i.  e. , 
the  classical  critical  speed  calculation). 


E)  The  bearing  reaction  F  lbs. 

3)  The  rotor  stiffness  k  ^0/in.  Since  the  results  are  valid  only  for  a  two 
bearing  rotor  and  for  rotor  speeds  below  the  second  critical  speed, 
the  rotor  stiffness  ia  calculated  from  k  =  M  <  c ol  where  M 

is  the  vibratory  mass  of  the  rotor,  M  is  somewhat  smaller  than  the 
actual  rotor  mass  and  may  be  estimated  by  methods  as  shown  in 
"Vibration  Problems  in  Engineering",  by  Timoshenko,,  Chapter  1, 
Article  4. 


4)  The  journal  radius  R,  in. 


For  later  use,  we  shall  define  the  bearing  parameter  X  as; 


X 


lbs  •  sec 


( jjt,  ;  viscosity,  —  -  P:  bearing  radius,  in  »  L;  effective  bearing 

length,  in  -  C;  radial  bearing  clearance,  in.) 
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In  addition  wo  shall  define: 


the  dimensionless  hearing  reaction: 


F 


the  dimensionless  rotor  stiffness 


the  dimensionless  speed  ratio 


(u»i  operating  rotor  speed, 


rad 


/sec) 


Once  these  three  parameters  are  known  the  transmitted  force  can  be  found 
directly.  Thus,  the  problem  is  to  choose  the  bearing  dimensions  in  such  a 
way  that  the  value  of  the  above  parameters  minimise  the  transmitted  fores, 


J± 

Selection  of  bearing  type  and  of  D-  ratio 


The  selection  will  be  done  on  a  basis  of  comparison,  For  this  purpose 
it  is  necessary  first  to  estimate  a  bearing  clearance  C  and  an  oil  viscosity  /w  . 
Since  the  journal  radius  R  is  known  It  is  then  possible  to  compute  A  for 
■jjy  a  -Ij-  and  _jjji.nl,  In  addition  and  are  obtained.  For 

a  given  operating  speed  w  ra^/sec  the  transmitted  force  can  be  found  as 
follows: 


id, 

a)  calculate  the  speed  ratio  wt  . 

b)  calculate  and  enter  figure  £  -10  to  get  the  corresponding  value 

of  eccentricity  ratio  g  , 

c)  enter  figure  11  with  ^  to  find  the  equivalent  speed  ratio  (^j  . 

d)  enter  figure  28»<*3  with  fcfif  and  the  corresponding  value  of  6  to  find 

the  dimensionless  transmitted  force  for  all  desired  bearing  types  and 


D 


ratio.  If  the  curves  are  spaced  too  far  for  linear  interpolation  it  is 


necessary  to  make  a  cross-plot.  Multiply  the  result  by  e  :  distance 

between  shaft  center  and  center  of  gravity  of  rotor  mass,  inches,  e  may 


be  calculated  from  the  equation:  c  = 

to  obtain  the  transmitted  force  P  lbs, 
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rotor  unbalance,  lbs  :  in 


‘/rotor  weight,  lbs.) 


This  procedure  can  be  repeated  for  a  number  of  operating  speeds  to 

cover  the  entire  operating  speed  range.  By  plotting  the  curves  of  transmitted 

Xj 

force  versus  rotor  speed)  the  effect  of  bearing  type  and  of  -jj~  -ratio  is 
readily  seen  and  a  selection  on  the  basis  of  minimum  transmitted  force  can 
be  made.  An  example  of  the  results  obtained  by  the  outlined  method  is 
given  in  figures  44-45. 

Selection  of  bearing  clearance 

In  principle,  the  selection  of  clearance  is  done  by  the  same  procedure 
as  above.  Thus  the  goal  ia  to  obtain  curves  of  transmitted  force  versus  rotor 
speed  for  various  values  of  the  bearing  clearance  and  from  that  select  the 
actual  clearance  value  on  the  basis  of  minimum  transmitted  force. 

As  before,  the  rotor  is  assumed  known.  In  addition,  an  oil  viscosity  must 
be  chosen,  Since  the  journal  radius  R  is  given,  it  is  then  possible  to  compute 
A  and  consequently  ^  and  for  JL »  Jp»  and  »  1  and 

for  various  values  of  the  clearance  C. 

As  it  would  be  time  consuming  to  cover  the  complete  speed  range,  it 
should  be  sufficient  to  base  the  comparison  on  rather  few  points.  This  is 
most  easily  done  in  the  following  way! 

a)  enter  figure  8-10  with  &«,2,  ,5  and  .7  (except  for  the  elliptical  bearing 

with  ellipticity  m  =,  5  where  the  values  are  &  a.  15,  .  3  and  .  5)  to  find 

r 

r<zx  • 

b)  divide  the  result  into  to  get  and  enter  figure  11  with  ^  to 

obtain  the  equivalent  speed  ratio  ($.) 

c)  enter  figures  28-43  with  the  values  of  ($■)  to  the  intersection  with  the 

corresponding  e -curve.  Read  off  the  value  of  and  multiply 

by  to  find  the  actual  transmitted  force  P  lbs. 

Jb 

Thus,  for  each  clearance  value  and  for  each  bearing  typo  and  D  -ratio, 

three  points  are  readily  obtained  on  the  curve  of  transmitted  force  versus 
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rotor  speed.  A  fourth  point  is  P»0  for  *0, 

Although  four  point®  are  not  sufficient  te  define  the  graph  of  P  versus 
speed  with  any  high  degree  of  accuracy,  it  may  at  least  be  enough  to  serve 
as  a  basis  for  comparison  and  a  subsequent  selection  of  bearing  clearance, 

The  effect  of  oil  viscosity  is  treated  in  the  same  way, 
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DYNAMIC  OPERATION 


Attenuation  of  structure  born*  noise  through  bearings  is  only  one  dynamic 
characteristic  o£  bearings,  Thera  are  several  other  examples  of  dynamic 
operations,  a.g.i 

a, )  Transient  conditions  during  starting  or  shutting  down. 

b. )  Orbiting  of  the  journal  under  a  constant  vibrating  load, 

a.)  Motions  of  the  journal  center  under  oscillating  loads. 

The  locus  of  the  journal  may  be  in  a  closed  path  of  fixed  amplitude  or  it 
may  increase,  In  the  latter  case  the  system  may  be  unstable  causing  the 
journal  to  rub  the  bearing.  Fig,  50  shows  savsral  cases  of  dynamic  behavior. 

The  previous  section  which  deals  with  tha  analysis  of  structure  borne  noise 
through  bearings  indicates  that  for  most  effective  noise  attenuation  tha  journal 
should  operate  at  low  eccentricity  ratio.  On  tha  other  hand  It  has  alao  been 
shown  (Rtf.  1)  that  rotors  operating  at  low  eccentricity  ratio  are  susceptible 
to  instability  at  relatively  low  speed,  In  fact  it  has  been  proved  theoretically 
and  experimentally  (Ref.  2  )  that  the  threshold  of  instability  for  vertical 
rotor  in  plain  cylindrical  journal  bearings  is  zero  speed.  These  two  conditions 
are,  therefor*,  somewhat  incompatible.  Since  the  system  must  be  stable  it  is, 
therefore,  necessary  to  optimise  noise  attenuation  without  sacrificing  stability, 

To  make  the  report  more  complete  a  section  dealing  with  Stability  and 
Balancing  is  also  diaeuased  briefly  in  this  report,  Several  definitions  are 
given  so  as  to  clarify  the  usage  of  some  of  the  terms  used. 


DYNAMIC  OPERATION 


* 


under  this  condition,  the  journal  center  moves  relative  to  the  bearing 
center  and  tha  local  fluid  film  pressures  vary  with  time.  (Ref,  3) 

la&fcsfeiAia 

This  1»  a  dynamic  condition  In  which  the  journal  center  moves  away  from 
the  bearing  canter,  until  breakdown  of  the  film  occur?  and  there  is  physical 
contact  between  the  journal  and  the  bearing,  Another  condition  that  mey  be 
defined  as  unstable  Is  that  In  which  the  Journal  canter  whirls  along  a  random 
locus. 

Threshold  of  Instability 

Corresponds  to  frequency  at  which  Instability  is  Initiated. 

MisosaObie 

This  Is  a  resonant  vibration  of  a  journal  In  a  fluid-film  bearing  which, 
for  low  eccentricity  ratio,  sets  in  at  approximately  twice  the  actual  first 
system  critical  and  persists  at  higher  speeds  with  frequency  of  vibration 
approximately  equal  to  the  first  system  critical  rogardleas  of  running  speed, 
The  motion  of  the  shaft  center  is  in  the  same  direction  as  shaft  rotation. 
Resonant  whip  Is  a  self- supported  vibration,  aa  is  half- frequency  whirl,  In 
the  case  of  resonant  whip,  the  vibration  Is  supported  by  the  fluid  film  action, 
while  the  frequency  la  controlled  by  the  system  critical  speed. 

Critical  Speed 

Critical  speed  Is  the  rotating  speed  of  a  system  which  corresponds  to 
resonance  frequency  of  the  system.  The  system's  critical  speeds  Include 
rigid  body  as  well  aa  bending  or  torsional  critical  speeds.  (In  this  text 
when  we  refer  to  first  critical  we  mean  bending  body  critical.) 
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Synchronous  Whirl 

This  is  a  whirling  orbital  motion  of  the  journal,  at  a  frequency  equal  to 
the  rotational  frequency.  The  motion  of  the  journal  is  in  the  same  direction 
as  the  direction  of  the  rotating  member. 

An  example  of  the  synchronous  whirl  is  the  case  of  unbalanced  rotating 
load.  (In  the  case  of  vertical  rotor  in  plain  cylindrical  journal  bearings, 
the  whirling  locus  is  a  circle!  in  the  case  of  horizontal  machine  with  plain 
cylindrical  Journal  bearings  the  whirling  locus  la  an  ellipse.  See  Fig.  SO,) 


STABILITY 

Using  the  coordinate  system  of  Fig.  51  the  dynamic  equations  m«y  be 
represented  by 

a.  Jk 


. .  OhL  + 

eoG 

eoC 

Se 

C  r 

“>£J 

1  +  Z£rt 

to  9f 

r  n  . 

«oC 

■ +  7tr~ " 

0 

Se 

C 

t  tv 

TT  a 


Sf 


f  (*  +  £) 


M  ...  .  ' \ 


kx 

2 


ky 

2 


(31) 


Here  the  functions  £„,  £.  and  their  derivatives  with  reapact  to  e,  e’  are 

4  C 

all  evaluated  at  the  equilibrium  eccentricity  ratio  c  .  and  for  c1  «  0. 

de  ■  ,  da  ■  -“-s 

C  eoC. 


\  ULR  R 
k  »■  sz^r  q 


,  M  is  rotor  mass, 


«  Id  - -  .  ,  £ 

-  “  %r  "  o5  ’  d0i  "  Tv 

0 

The  differential  Eqs.  (31)  are  linear  in  the  variables  £,  q;  x*  y,  and 
their  solutions  contain  time  as  an  exponential,  These  may  bo  expressed  in 
dimensionless  form  as 
„VT 


e ' '  where  t  »  «  fc 
o 


01  K 
0 


(32) 
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1 


Here  cd&  Is  the  critical  speed  of  the  simply  supported  shaft-rotor  system 
whose  mass  Is  M  and  whose  stiffness  Is  k, 

From  the  right-hand  pair  of  Eqis.(3i)f  there  results 


M  (vw0)' 


k  +  M  (VCD.)' 


s*  M  <«“„>'  „ 

£  y  ■  -  77 

k  +  M  (vcd  ) 2  1 


(33) 


Equation  (33)  now  leads  to  the  determinants!  equation 
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'  C  0  r/\  06  o  oc  2  \|k  +  M  (w >o)z 


If  w#  now  introduce  the  dimensionless  ratio 


■  0 


(34) 


(D_ 

CD. 


(35) 


where  cd  Is  the  angular  speed  at  the  threshold  of  instability,  we  get 

,Z  ■  •  c3£  &£. 
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(37) 


By  writing 


r- 


Av 


t  0.  \/ 


(38) 
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Eq, (3fi)  becomes 


2 

Is  not,  in  general,  equal  to  aero,  so  that  the  factor  can  be  divided 
out  of  Eq.  (39), 

It  was  assumed  in  the  derivation  of  Eq, (34)  that  the  solutions  of  the 

VT 

equations  of  motion  were  of  the  form  e  ,  where  V  is  a  complex  number .  If 
the  system  is  dynamically  stable,  the  real  part  of  the  complex  number  V  is 
negative.  Conversely,  if  the  system  io  dynamically  unstable,  the  rati  part 
of  v  is  positive.  Thue,  at  the  threshold  of  instability,  v  will  be  a  pure 
imaginary  number. 

We  now  solve  Eq,(39)  for  the  condition  where  v  is  wholly  imaginary  in 
order  to  obtain  the  value  of  u>  at  the  onset  of  instability, 

Considering  first  the  imaginary  part  of  Eq. (39),  we  have 
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If  8  »  0,  we  obtain  a  trivial  solution;  for  s  +  0,  we  have 
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Hast  considering  the  real  part  of  Eq.  (39),  m  have 
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From  Eq.  (38)  we  have 

fv2  -  Av2  +  £  -  0 
Once  again,  for  a  r  0,  we  can  write 
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The  speed  «  at  which  instability  starts  to  occur  is  now  defined,  since 


«)  «  SO)  . 

o 

The  above  defined  speed  at  which  instability  sets  in  is,  in  general, 
different  from  the  critical  speed  of  the  shaft-rotor-bearing  system.  For 
a  symmetrical,  two-bearing  system  the  critical  speed  may  be  calculated  as 
follows ! 
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The  critical  spaed  of  the  system  is  then 


where  subscript  v  refers  to  rsdiel  stiffness. 

The  dimensionless  number  A  [defined  in  Bq .  (37)]  ,  is  e  function  of 
beering  geometry,  shaft  stiffness,  and  fluid  viscosity.  Calculations  for 
the  threshold  of  instability  in  which  A  was  varied  from  0.1  to  100  for  0.1  <£  c  < ^ 
0.8  and  L/D  ■■  0.5  and  1  were  performed.  The  values  of  f  ,  £.,  df/de,  Sf./Sc, 

Sf  /Sc*  and  Sf„/Sc'  were  obtained  from  the  solution  of  the  dimensionless 
Reynolds  equation.  By  Introducing  these  values  into  Iqa,Q7),  (40),  (42),  and 
(43),  ws  obtain  the  results  of  Table  6.  The  results  indicate  that,  while  for 
low  eccentricity  ratios  instability  sets  in  at  approximately  twioe  the  critical 
speed,  this  number  increases  with  an  Increase  in  eccentricity  ratio,  Thus,  the 
onset  of  instability  for  eccentricity  ratios  of  0,8  is  about  four  times  the 
critical  speed,  This  conclusion  agrees  with  observations  which  show  that 
stability  increases  with  an  increase  in  eccentricity  ratio  and  also  that 
instability  may  occur  even  at  high  eccentricity  ratios. 

•mmummm 

The  number  (l/i)(v/s)  shown  in  Table  ^  (where  i  «  v  '1  )  represents  the 
ratio  of  the  frequency  of  the  oscillation  of  the  shaft  center  to  the  running 
frequency  of  the  shaft,  calculated  at  the  onset  of  Instability,  Note  that  this 
ratio  is  always  below  0.5  and  is  independent  of  the  magnitude  of  A. 
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BALANCING 

A  rotor  is  said  to  be  balanced  perfectly  when  it  rotates  in  free  space 
about  one  of  its  principal  axes  of  inertia,  which  would  be  an  axis  of  symmetry 
if  such  exists,  with  no- wobble.  H  circular  journals  are  now  constructed 
concentric  with  this  axis,  tey  definition,  these  journals  would  also  rotate 
with  no  wobble  and  may  be  enclosed  in  bearing  housings  with  no  rotating  forces, 

From  this  it  follows  that  journals  which  are  not  concentric  with  a 
principal  axis  of  inertia  tend  to  wobble  and  for  a  rotor  to  revolve  about  «n 
axis  defined  by  such  journals  it  must  be  driven  by  a  set  of  applied  forces 
which  rotate  with  the  rotor  If  the  center  of  gravity  of  the  rotor  is  not  on 
the  lino  of  bearing  center  there  must  be  a  net  force  which  furnishes  the 
centripetal  acceleration  of  the  center  of  gravity.  This  force  is  equal  to  the 
maea  of  the  rotor  times  the  centripetal  acceleration  of  the  G,G.  The  C.G.  may 
be  moved  to  the  axis  of  bearing  centar  by  suitable  weight  or  weights.  Such 
correction  is  commonly  called  static  balancing. 

Once  the  C.G.  is  moved  to  the  line  of  bearing  centers  ft  rotating  couple 
must  be  applied  to  keep  the  rotor  rotating  about  an  axis  at  on  angle  to  that 
of  «  principal  axis  of  inertia.  This  couple  may  be  rxorted  by  the  centrifugal 
action  of  two  equal  and  opposite  weights  in  any  two  arbitrary  planes.  This 
correction  plua  that  of  the  G.G.  correction  is  known  as  dynamic  balance  and 
may  only  be  determined  on  a  rotating  rotor-bearing  system. 

The  most  fundamental,  but  not  unusual,  description  of  unbalance  appears 
to  be  that  of  the  deviation  of  the  line  of  bearing  centers  from  the  principal 
axis  of  inertia.  It  is  important  to  note  that  the  specifications  of  unbalance 
depends  on  the  location  of  the  line  of  bearing  centers.  The  two  simplest 
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descriptions  of  unbalance  would  appearto  be  the  displacement  of  the  C.G,  and 
the  angle  of  the  inertia  axis  from  the  bearing  axis  or  the  displacement  of  the 
geometric  canters  of  the  journals  from  the  bearing  axis.  (Refer  to  Fig.  SI) 
Under  dynamical  conditions  either  e  or  a  or  both  may  not  be  constant.  The 
conditions  of  dynamical  equilibrium  gives  rise  to  the  following  equations. 
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where 


cb  »  is  the  angular  acceleration  of  the  rotor  (generally  very  small) 
8  ■  is  distance  between  rotor  geometric  center  and  mass  center 


These  are  equations  of  motion  with  e  and  a  as  the  two  degrees  of  freedom. 
Once  the  fluid  film  forces  are  known  it  is  possible  to  solve  these  equations, 
The  fluid  film  forces  are  to  be  found  by  integrating  the  pressure  over  the 
projected  journal  surfaces  normal  and  parallel  to  the  plane  of  maximum  film 
thickness  respectively. 
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The  fluid  film  pressure  satisfies  the  generalized  Reynolds  aquation 
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Thus  it  can  b«  shown  that  the  radial  and  tangential  forces  are  a  function  of 

■  ^ (S)  (m •  ii <«»  tuHtor 1 1/°)  (5W 

'  »t  -  Jr  jlf  c»  -  «a>  »t  <«,  ,  wd>  (57b> 

«  * 

Thus  for  a  specified  bearing  geometry  and  speed  of  rotation  once  c,  c  and  & 
are  measured  the  fluid  film  forces  can  be  readily  calculated.  They  may  be 
constant  per  cycle  e  »  0  (e.g.,  vertical  rotor  in  plain  cylindrical  journal 
bearings)  or  they  may  vary  from  point  to  point  along  the  journal  locus  this 
case  corresponds  to  the  horizontal  rotor  with  gravity  and  unbalance  load.  Once 
these  forces  are  established  the  magnitude  and  phase  angle  of  the  unbalance 
force  can  be  calculated  from  the  dynamic  equations  (54a, b) .  Thl*  permits  balancing 
of  tha  rotor  without  *  trial  and  error  procedure. 

* 

The  measurements  of  displacement  (a)  and  velocity  (c,a)  of  journal  center 

»  can  be  obtained  by  use  of  two  capacitive  or  inductive  probes  located  at  90°  to 

each  other  within  the  bearing  bore.  They  would  then  measure  the  motion  of  the 
journal  center  with  respect  to  the  bearing  center  as  a  function  of  time.  These 
capacitive*  or  inductive  probes  can  also  serve  as  monitoring  devices  to  determine 
bearing  performance  in  service.  Thus  they  serve  a  dual  function  of  providing 
the  necessary  measurements  for  balancing  and  monitoring  bearing  performance. 

Figs.  48  and  49  show  the  locus  of  the  shaft  center  that  such  pickups  would  aoo. 

Without  the  use  of  such  instrumentation  within  the  bearing  there  Is  still 
another  method  which  may  be  employed  for  balancing.  The  analysis  presented  in 
this  report  indicates  that  there  is  a  relation  between  the  driving  force  and 

* 

the  transmitted  force.  The  difference  being  the  attenuation. 
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Since  She  attenuation  has  been  theoretically  calculated  by  meaauramant  of 
force  transmitted  It  Is  possible  to  calculate  the  driving  force  and  the  phase 
angle.  Once  this  is  done  it  is  possible  to  balance  the  rotor  without  resorting 
to  trial  and  error  procedure, 

This  analysis  indicates  that  by  either  measurements  of  journal  locus  or 
transmitted  force  it  ia  possible  to  determine  the  magnitude  and  phase  angle  of 
the  unbalance.  Corrections  sen  then  be  incorporated  to  balance  the  system, 

To  illustrate  this  point,  theoretical  and  experimental  analysis  have  been 
carried  out  on  synchronous  whirl  (a.g.,  whirl  produced  at  running  speed, 
unbalane*'  load)  with  circular  orbit  about  the  bearing  center  (eug,,  vertical 
rotor),  (Raf,  4  )  The  comparison  between  the  measured  unbalance  force  and 
phase  angle  and  calculated  values  appear  to  be  very  good,  similar  comparisons 
will  b@  carried  out  on  horizontal  rotor, 


coKcmsi-OKrg 


1,  Attenuation  of  structure  borne  noise  should  be  done  as  cioB.s-to  the  source 
of  noiee  aa  nosaihle,  Fluid  film  bearings  provide  the  source  of 
attenuation. 


Qjhn a 


t.  The  report  provides  means  for  calculating  noise  attenuation  of 
rotors  supported  by  bearings  of  three  geometries.  The  results  are 
presented  in  dimensionless  form  ao  aa  to  be  applicable  to  a  large 
range  of  geometrically  similar  configurations.  The  report  also  pro¬ 
vides  means  for  calculating  system  critical  speeds  and  the  static 
load  carrying  capacity, 


3.  The  report  shows  that  fluid  film  bearings  can  provide  considerable 
amount  of  viscous  damping  which  absorb  the  vibrational  energy  and  in 
this  way  attenuate  the  force  transmitted  to  the  structure.  Bearing 
geometry  plays  a  major  role  on  the  level  of  attenuation  as  exemplified 
by  the  differences  of  the  three  bearings  studied,  see  figures  44»45. 

4.  Since  the  fluid  film  is  an  elastic  media,  accurate  predictions  of  the 
system  critical  speeds  must  include  the  elasticity  and  the  damping 
of  the  Journal  bearings,  see  figures  46 ~47 . 


5.  Indications  are  that  for  maximum  attenuation  one  should  operate  at 
low  eccentricity  ratio,  low  stiffness  and  high  damping.  However, 
this  implies  low  critical  speed  anHljenHeney  for  resonant  whip. 
Therefore,  one  must  optimize  the  design  to  ensure  stability  and  at 
the  same  time  get  maximum  attenuation. 


6.  The  dynamic  response  of  a  rotor  can  bo  measured  bv  capacitor  or  ^  kv 

-inductive  pick-upg .  These  can  serve  also  as  monitory  devices  for 
controlling  bearing  pe rformance.  Obtaining  the  dynamic  response  and 
knowing  the  bearing  characteristics,  it  is  possible  to  determine  the 
magnitude  and  phase  angle  of  tho  unbalance  load  thus  eliminating  trial 
and  error  in  balancing. 


1 
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j^CO  wivi  SNB  A  T  IONS 

1,  Other  geometries  aueh  as  pivoted  shoes  and  elastically  supported 
members  in  which  stiffnsas  and  damping  can  be  controlled  and  varied 
should  be  investigated.  These  studies  should  be  both  theoretical  and 
experimental  for  these  geometries  offer  considerable  potential  in  noise 
attenuation  and  at  the  same  time  have  a  tendency  of  being  stable. 

Z,  Further  studies  should  be  continued  to  determine  rotor  dynamics 
around  the  second  bending  critical  and  higher. 

3.  Computational  techniques  should  be  set  up  which  would  include, 
without  transformations, the  bearing  cross -coupling  coefficients. 

4.  The  application  of  capacitive  or  inductive  pick-ups  as  permanent 
measuring  devices  in  bearings  should  bo  investigated.  Such  measure¬ 
ments  provide  means  for  continually  chocking  rotor  balance  and  for 
monitoring  bearing  performance. 
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0 
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0 
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,17457 
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.7 
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18.44.3 
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17,607 
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37,404 
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0 

68,197 
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37,231 

0 
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.90262 
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.13 

B0 
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0 

1,6499 
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,91698 
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,(19049 
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0 
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DutYHimg  reaffirm  m  a^iristim.  <o t  clocity  in  y-dimiimt 
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Dkm^ifii  Goffitiant  m  y^iraeiiHK  for  valoiiy  In  y«dirasilon 
Snaring  diarnawr 
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Soaring  roetnm 
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tUkfina  fore*  m  miu.l  and  h*ria«Mal  airs*  lien 
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Oimgniionlo*  oil  tilm  tliuxnoi 
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H«<9*  HI  Hum* 

Kffoiiva  boring  langili 
Vibratory  rotor  mail 

Roaring  tlliimeiiy,  in  (igum  I 
Rotor  »i»«(f 
Oil  film  pr«MUH 

Pure*  iranimiiitd  by  b*»rlog  in  rarlifiil  ami  Noriiontal  dlraclinn 
tarri  irantmitiad  for  ngnl  rotor  lopporii 
n»*ring  radiu* 

Clrci»mf*r*niial  ami  axial  conriltnats*  for  ml  film 
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Qiv«n  by  m,  (III,  pig*  n 
Raaring  octanirimiy  ratio 
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I’olar  roordlnaln  for  oil  film 
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